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We show that Brownian motion of a one-dimensional domain wall in a large but finite system yields 
a u -3 / 2 power spectrum. This is successfully applied to the totally asymmetric simple exclusion 
process (TASEP) with open boundaries. An excellent agreement between our theory and numerical 
results is obtained in a frequency range where the domain wall motion dominates and discreteness 
of the system is not effective. 
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The asymmetric simple exclusion process (ASEP) is 
one of the simplest models of collective transport of 
particles. Q, 0> Ej It can be interpreted as a simplified 
model of traffic flow 4, 5] or as a model of granular flow. |(J 
Since some exact solutions of ASEP can be obtained un- 
der given boundary conditions, we may capture essential 
feature of more complicated systems such as traffic flows 
and granular flows. 

In the study of the traffic flow or granular flow, the 
power spectrum I(ui) is used to investigate properties of 
local density fluctuations of vehicles or particles and it 
commonly shows a power-law form I(uj) ~ uj~ a . For ex- 
ample, in real observation of granular flow, it is confirmed 
that a = 4/300,l3. On the other hand, the value of the 
exponent a has not yet fixed for traffic flows. 0, ^} It is 
worthwhile mentioning that the structure factor (Fourie 
transform of the power spectrum) in the Rouse model, 
which is the simplest model of polymer dynamics, obeys 
exp[— (t/r) 1 / 2 ] with the relaxation time r 11]. Namely, 
the Brownian motion of structured material produces 
I{lj)~uj- 3 / 2 . 

In this paper, we numerically calculate the power spec- 
trum for the totally asymmetric simple exclusion process 
(TASEP) which is the simplest ASEP and find that this 
model also exhibits power-law I {to) ~ w~ 3 / 2 . This be- 
havior can be explained by a domain wall theorv p^.ll3| . 

The TASEP is a continuous-time stochastic process 
on the one-dimensional lattice which is defined as fol- 
lows. Consider the set of lattice points —L/2, —L/2 + 
1, . . . , L/2. Each site i is either occupied by a particle 
(tj = 1) or empty (t; = 0). During an infinitesimal time 
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interval dt, each particle hops to the right adjacent site 
with probability dt if the destination site is empty. More- 
over, in case of open boundaries, a particle is injected to 
site —L/2 with probability adt if the site is empty, and 
a particle is extracted from site L/2 with probability (3dt 
if the site is occupied. 

The steady state of this model is exactly obtained for 
any parameter values a and (3 an d sy stem size L by the 
method of matrix-product ansatz hjL which can also be 
applied to the more general ASEP[lja]. Moreover, in the 
steady state, the mean occupation at site i, p = (t;), 
and the current J = (tj(1 — tj+i)) can be calculated in 
the large L limit. As the result, it turns out that there 
exist three phases depending upon values of a and (3 as 
follows. When a < i and a < (3, the mean occupation is 
p = a and the current is J — a(l — a). This is called the 
low density phase. When (3 < | and (3 < a, the mean 
occupation is p = 1 — f3 and the current is J = (3(1 — (3). 
This is called the high density phase. When a > | and 
[3 > i, the mean occupation is p = k and the current is 
J = j. This is called the maximum current phase. The 
phase diagram is illustrated in Fig. 1. 

The thick line in Fig. 1, a = /3 < |, is the coexis- 
tence line, where we focus our attention in the following. 
In this case, there occurs a domain wall between a left- 
hand region of the low density phase and a right-hand 
region of the high density phase. It is known that the 
time-averaged density profile is a linear function of the 
position. This means that the domain wall can be any- 
where with equal probability as a result of its random 
walk. This kind of kink motion can also be obtained 
from a decouple approximation of ASEP where the equa- 
tion of density profile is reduced to the inviscid Burgers 
equation. It should be noted that the inviscid Burgers 
equation is a model equation of traffic flow0,|jj. 

The average motion and fluctuations of a domain wall 
is obtained as follows fl^. Il3|. Assume that the system 
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1: The phase diagram of the TASEP with open bound- 
aries. Regions A, B and C mean the low-density phase, the 
high density phase and the maximum current phase, respec- 
tively. 



consists of a left domain of particle density p_ and flux 
j_ and a right domain of particle density p + and flux j + . 
Then the drift velocity of the domain wall is derived via 
the equation of continuity as 



where u n 
trum I (id) is computed by 



Then the power spec- 



I(u n ) = T(\c/ )n \ 2 ), 



(5) 



where the brackets mean averaging over 256 samples. 

Figure 2 shows the simulation result for y = 0, which 
lies in the center of the system. It is clearly observed 
that the power spectrum is propotional to w -3 / 2 in the 
frequency range 5 x 10" 5 ~ 2 x 10' 3 . In the follow- 
ing, we show that this power law behavior is induced by 
Brownian motion of the domain wall. 
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This is further interpreted as the result of random walk 
of the domain wall with hopping rates 



D + = 



J+ 



D_ = 
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for a move to the right (D+) and the left (D-). In partic- 
ular, on the coexistence line, the drift velocity vanishes 
and the diffusion constant for the domain wall motion is 
given by 



D 



a(l — a) 
1 -2a ' 



(3) 



We carried out Monte Carlo simulations of the TASEP 
with open boundary conditions of system size L = 200. 
In our simulations, one Monte Carlo step corresponds to 
time interval 0.1. Parameter values are chosen to be a = 
(3 = 0.2, which is on the coexistence line. After discard- 
ing transient parts before a sufficiently large time T , se- 
quences of site values {x t — T y (T + t)\t = 0, 1, . . . ,T— 1} 
at a fixed position y is recorded. The length of the se- 
quence is T = 2 20 . Fourier components of the series is 
denoted by <f> n . That is, 



FIG. 2: Power spectrum fo time sequence {ro(t)} in the 
TASEP with open boundaries. 

Now we consider a Brownian motion of a domain wall 
in an idealized form. Let a domain wall be described by 
the step function 0(x) as 



r(x,t) = 6(x-X(t)), 



(6) 



where X(t) is a Brownian motion in [—L/2,L/2], The 
end points — L/2 and L/2 are assumed to be the reflecting 
boundaries. Then, the autocorrelation function of r(y, t) 
is given by 

rv rv 

(T(y,t)T{y,0)} = dx dxP st (xo)P(x,t\x Q ,0) 
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where P s t(x) is the stationary probability density for 
X(t) = x and P{x, t\xo, 0) is the transition probability 
density. In other words, if D denotes the diffusion con- 
stant for X(t), P s t{x) is the stationary solution of the 
diffusion equation 



dP _ D d 2 P 
dt dx 2 



(8) 



with the boundary condition 
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and P(x, t\xo, 0) is the solution with initial condition 
P(x,0\xq,0) = S(x — Xq). The equation is easily solved 
and we find 



Pst(x) = 



1 



and 
P(x,t\x ,0) 

where A„ = 
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with n = 1, 2,3, ... . Then, the integral 



in the right side of Eq. (JJJ is carried out to produce 

(r(y,t)r(y,0)) = 
(L + 2y) 2 °° -- DA -* 
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This result can be generalized to negative t by replacing 
t by \t\ in the above expressions. 

Because the power spectrum is equal to the Fourier 
transform of the autocorrelation function, we have 
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If the oscillatory part including factor (— 1)" _1 can be ig- 
nored and the sum over n can be replaced by an integral, 
we can evaluate it as 
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Thus we have arrived at the w -3 / 2 law. It is remarkable 
that this result does not depend on position y. Actually, 
as we will see later, I(to) shows deviations from the above 
result in small u) if \y\ is large. This is because the os- 
cillatory part cannot be ignored in case of large \y\ and 
small w even if system size L is large. 

Similar calculations can be done if the space is discrete. 
Assume that a domain wall performs a continuous-time 
random walk on the one-dimensional lattice between two 
reflecting walls located at — i and 4 . Denoting hopping 
rate by D, we obtain 




(L + l) 



,(15) 



where [y] denotes the largest lattice point that does not 
exceed y. As in the previous case, if the oscillatory part 
can be ignored, replacement of the sum in the right side 
by the integral yields for uj ^ 
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In case u> <C 4D, this returns to the previous result 



(16) 




FIG. 3: Power spectra of the time sequences of site values 
{r(y,t) = % - W(t))\t = 0, 1, 2, . . . } for y = and y = 0.3L, 
where W(t) means a random walk with hopping rate D = 0.1. 
The "theory (y — 0.3L)" means the right-hand side of Eq. 
115H and the "asymptotic form" means the right-hand side of 
Eq. (fTo) . 

To confirm the above results, we carried out numerical 
simulations. In the simulation, we deal with a system 
of length L = 200 and let a domain wall carry out a 
random walk W(t) with hopping rate D = 0.1. Figure 3 
shows numerically obtained power spectra in case y = 
and in case y = 0.3 L. In case y = 0, the agreement 
between the theory and the numerical result is excellent. 
Deviations from the power law are observed only in the 
very small frequency range ui < ir 2 D/L 2 ~ 2.5 x 10~ 5 , 
where discreteness of the spectrum {A„} is eminent. In 
case y = 0.3L, deviations extend to larger frequencies. 
As mentioned earlier, this is because the oscillatory part 
cannot be ignored if \y\ is large. In any case, however, 
the power law behavior is still clearly observed. 

Now we return to the TASEP on the coexistence line. 
Within the domain wall theory, t v (t) is approximated by 



r v (t) 



(1 



■0)6(y - X(t)), 



(17) 



where X(t) denotes the position of the domain wall 
at time t. Then, the prefactor of w~ 3 / 2 becomes 
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y "21^ 2 "' ) • Actually, the straight line in Fig. 2 
is depicted with this prefactor. The excellent agreement 
between the theory and the numarical result confirms 
that the power-law behavior in the power spectrum is 
induced by the Brownian motion of a domain wall. 
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FIG. 4: Power spectra multiplied by system size for the 
TASEP of various system sizes: L = 100, 200 and 500. 

The prefactor of the power law includes inverse of sys- 
tem size L . This means that the power law in the 
power spectrum is a finite-size effect and is observed only 



in large but finite systems. Figure 4 shows numerical re- 
sults with various system sizes. As system size L varies, 
the power-law part of the power spectrum is scaled 
by L _1 as expected, while higher frequency part is not. 
In fact, this part is independent of system size L. This 
means that the higher frequency part represents local 
fluctuations. 

In summary, we have derived that Brownian motion 
of a domain wall yields an oj~ 3 / 2 power spectrum. This 
mechanism is successfully applied to explain the power 
law in the power spectrum for local density fluctuations 
in the TASEP. 

It should be noted that the ASEP does not admit more 
than one domain wall. In real traffic or granular systems, 
however, many domains can occur and interactions be- 
tween them may be important. It is a future problem 
to investigate power spectra in such more complicated 
systems. 

It is remarkable that the domain wall theory is quan- 
titatively valid for TASEP. This success may suggest the 
flutcuation of kink position is characterized by a simple 
diffusion process. On the other hand, it is known that the 
current fluctuation is not Gaussian[jj|. We will have to 
clarify this puzzled situation in the fluctuations of ASEP 
in the future. 
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